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Abstract 

The Broad Emission Lines (BELs) in spectra of type 1 Active Galactic 
Nuclei (AGN) can be very complex, indicating a complex Broad Line Region 
(BLR) geometry. According to the standard unification model one can expect 
an accretion disk around a supermassive black hole in all AGN. Therefore, 
a disk geometry is expected in the BLR. However, a small fraction of BELs 
show double-peaked profiles which indicate the disk geometry. Here, we 
discuss a two-component model, assuming an emission from the accretion 
disk and one additional emission from surrounding region. We compared the 
modeled BELs with observed ones (mostly broad Ha and H/? profiles) finding 
that the model can well describe single-peaked and double-peaked observed 
broad line profiles. 
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1. Introduction 



The detailed structure of the innermost part of AGN is still an open 
question. It is widely accepted that the central engine consists of a super 
massive black hole fueled by an accretion disk. The origin of the broad 
emission lines, observed in most AGN, is explained by photoionization of 
dense gas surrounding the nuclei, but the structure and dynamics of the BLR 
remains unclear. Using Broad Emission Lines (BEL) shape investigations of 
the BLR, one can constrain the geometry of the material in the BLR (see 



Popovic 



20061 ). 



BELs of AGN show high diversity in their shapes and widths. In some 
rare cases, where the accretion disk clearly contributes to the BELs, we can 
invest igate the disk properties using a broad, double-peaked, low-ionization 



lines 



2003 



Chen fc Halpern. 



Eracleous et. al 



1989 



Chen et al. 



1989 



Eracleous fc Halpern 



1994 . 



20091 ) . The rotation of the material in the disk results 
in one blue-shifted and one redshifted peak, while the gravitational redshift 
produces a displacement of the center of the line and a distortion of the line 
profile. 

Moreover, a number of BELs show shoulders-like profile in t he wings of the 



line, which could be a fingerprint of the accretion d i sk em ission (IPopovic et al. 



2003 



Eracleous fc Halpern! . 



2003 



Strateva et al. 



20031 ) . Beside emission of 



the dis k (or disk-like region) or emission from 



a disk (IChen fc Halpern . 



1989 



Chen et al. 



spiral shock waves within 



19891 ). the following geometries 



may cause substructures in line profiles: i) emission from the oppositely- 
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directed sides of a bipolar outflow (jZheng et al. 



1990 



1991 



); ii) emission 



from a spherical system of clouds in rand omly inclined Kep l erian orbits il- 



luminated anisotropically from the center 



emission of the binary black hole system (jGaskelll . 



(IGoad fc Wanders 



1983 



1996|); 



19961 ): and hi) 



To explain the complex morphology of the observed BEL s hapes, different 



geom etrical models have been discussed (see in more details 



Sulentic et al. 



20001 ) . In some cases the BEL profiles can be explained only if two or more 



kinematically different emission regions are cons i dered (see e .g. 



1993 



2006 



Romano et al 



Ilic et al. 



2006 



996 



Popovic et al 



CoUin et al 



2006; 



2001 



2002, 



Hu et al 



2003 



Marziani et al. 



2004 



Bon et al 



20081 ). In particular, the 



existence of a Very Broad Line Region (VLBR) with random velocities at 
5000-6000 km/s within an Intermediate Line Re gion (ILR) has also been 



considered to explain the observed BEL profiles (ICorbin fc Boroson . 



Sulentic et al. 



2000 



Hu et al. 



1996 ; 



20081 ). 



Nevertheless, the majority of AGN with BELs have only single peaked 
lines, but this does not necessarily indicate that the contribution of the disk 
emission to the BELs profiles is negligible. 

In this paper we study the possibility that the flux of the hidden disk 
emission could be present in single peaked emission line profiles of AGN. 
With this goal, we applied a simple, two-component model, which assumes 
that observed BELs can be represented by composition of the emission of 
the accretion disk and surrounding non-disk isotropic reg ion. For the ac- 



cretio n disk emission, we have used the model proposed by I Chen &: Halpern 
( 119891 ). which assumes optically thick and geometrically thin disk, while for 
the additional emitting region which surrounds the disk, we assumed isotropic 
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velocity distribution of emitting clouds, which imply that the emission line 
profiles generated by this region can be described by a Gaussian function. In 
Section 2 of this paper we describe samples of galaxies where we analyzed the 
possible accretion disk contribution to the total flux of their broad emission 
lines. We review the methods of analyzes in Sections 3 and in Section 4 we 
give results. Finally, in Section 6 we outline our conclusions. 



2. Data sets 



First sample of 12 AGN spectra was observed with the 2.5 m INT on La 



Palma at 2002 (IPopovic et al. 



20041 ). For some lines with bad S/N ratio we 



used HST observations obtained with the Space Telescope Imaging Spectro- 



graph (STIS) on January 2000 



(NGC 3516 



data reduction can be find in (IBon et al. 



Details about observation and 



200J). 



2006; 


Popovic et al. 


. 2004 


The 


; e.g. 


Popovic et al.. 


2002. 


2003. 



With the idea to compare our model with real emition line profiles, with- 
out fitting, we used the sample of sp ectra containing 90 bro ad-line-emitting 



AGN, which have been collected by 



La Mura et al. 



data release of Sloan Digital Sky Survey (SDSS 



corrected for sky-emission, tel. 







(120071 ) from the third 



The spectra were already 



uric absorption, the Galactic extinction and 
redshift ( lLa Mura et all 120071 ). Since the interest was to investigate the 
broad emission line profiles, we subtracted the narrow components of Ha 
and the satellite [Nil] lines. The spectral reduction (including subtraction of 



1 http://www. sdss.org/dr3 
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stellar component) and the way to o 



details explained in 



La Mura et al 



jtain the broad line profiles are in more 



f|2007h . 



Previously cleaned broad Ha profile was normalized, converted from wave- 
length to velocity scale, using DIPSO software package. 

The sample of 90 AGN contained the spectra with different FWHM of 



broad H Q from about 100 km/s to 7000 



sm/s, with at le ast se veral represen- 



tatives in every 500 km/s lLa Mura et al. 



fsee Table 2 in 



20071). This covered 



many different types of broad single-peaked profiles. 



3. Methods of BEL profile analyzes 

We performed several tests to find a possible accretion disk contribution 
to the total line flux of a single peak broad emission lines. 



3.1. Gaussian fit 

We used a x 2 minimization routine to obtain the best fit. To limit the 



number of free parameters in t 



as given in iPopovic et al. 



re fit, we have set some a priori constraints 



(120031 ) 



In the fitting procedure, we looked for the minimal number of Gaussian 
components needed to fit the broad component. For the narrow blending 
lines we have been taking into account the intensity ratio of [OIII] lines ( 
obtained form the atomic values 1:3.03) and the fit of Fe II template for the 
Hp In the case of H a line, for the narrow [Nil] lines we assume that their 



2 http: / /www. starlink.rl.ac.uk 
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Figure 1: Decomposition of the Ha line of Mrk 841. The dashed lines represent the obser- 
vations, while the solid lines show the profile obtained by Gaussian decomposition. The 
Gaussian components are presented at the bottom. The narrow dashed lines correspond 
to the narrow Ha and [Nil] lines. 



intensity ratio is 1:2.96 (see 



Popovic et al. 



2003 . 



2004 



Bon et al. 



20061 ) 



Some examples are presented in Figures [T] and [2j 

The results of indicated existence of two kinematic regions, one very broad 
line region (VBLR), represented by the two shifted Gaussians and one in- 
termediate broad line region (ILR), represented with the central Gaussian 
component. It was noticed that both shifted components had similar widths, 
while the central component had similar values to the central components 
of the other lines (see Fig. 3), what led to idea to introduce the two com- 
ponent model, that consists of the accretion disk (which could describe the 
wings of the line), and surrounding region with isotropic velocities of clouds 
(representing the core of the line, that can be described by a Gaussian). 
These results encoura ged u s to test in the same manner the sample of 



12 AGN flPopovic et al. 



2004 ). After Gaussian decomposition the results 
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Figure 2: Decomposition of the H/3 line of 3c 273. The dashed lines represent the obser- 
vations while the solid lines show the profile obtained by Gaussian decomposition. The 
Gaussian components are presented at the bottom. The narrow dashed lines correspond 
to the narrow H/3, [Olll] lines and Fell shelf. 

showed again two distinct kinematic regions. The results of Gaussian analysis 
showed that Fe II lines had comparable widths to the central component 



see Fig. [3]), indicat ing that they origina te in the same r egion, w 



rat was 



previously shown by 



Marziani et al. 



(120031 ) and recently by 



Hu et al. 



(120081 1 . 



3.2. Fit with a two- component model 

For the disk emission, we adopted disk model developed by 



Chen fe Halpern 



( 119891 ). We fitted the profiles with composition of the disk and Gaussian 



profiles (see Figures 4-6). Starting value of the width of the Gaussian was 
considered to be the value of the width of the central Gaussian obtained 
from the fit with three broad Gaussians components, for the same object. 
All those lines were expected to have a common inclination, since they orig- 
inate in the same object, so this constrain was taken into account. 



Using these results as a starting point, we fitted with the two-component 
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Figure 3: The widths (w) as a function of the shifts of the broad Gaussian components 
obtained for the Ha and H/3 lines of our sa mple of AGNs (open circles) and the Fe II 



template (full triangles) (jPopovic et al 



2004). 



model the profile co mposed as averaged Ha and H/3 line of each AGN (see 
Popovic et allbooA . 

The results of the fit are presented in Figs 4-6, and corresponding param- 
eters of two component model are given in Tables 1-3 for 12 AGN. As one can 
see from Figs 4-6, the line profiles can be well fitted with a two-component 
model. Depending on the inclination, local broadening and the dimension of 
the accretion disk, a line from the disk can appear as single or double-peaked. 

The parameters obtained by fitting the line profiles are given in the Tables 
1-3 for 12 AGN. Here we presented several tests of fitting with two compo- 
nent model. We fixed the emissivity parameter as p=3 (Table 1) and also 
as p=2.5 (Table 2), with allowing all other parameters to be free. While 
searching for the maximal inclination of the disk, we obtained values of pa- 
rameters presented in Table 3. For the parameter that corespond to the 
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shift of the disk component, there are several possib 



motion of emitters in elliptical disks (jEracleous et al. 



c cx 



planations, like a 



19951 ). tidal perturba- 



tion of the disk around a supermassive black hole by a smaller companion 
(jEracleous et al.1 . 119951 ). or some binary black hole configurations in which 
BLR is dominated by the gravity of one black hole. The orbital motion of two 
black holes around their center of mass can result in radial velocity shift o 



the e ntire BLR (see for example 



Boroson fc Lauer 



2009 



Bogdanovic et al. 



20091 ). Moreover, one can expect different types of turbulenc e caused by 



magnetic field in the accretion disk (see e.g. 



Turner et al 



20031 ) that in the 



case of the line additionally contribute to the line broadening. The magnetic 



field could produce the effects of t 
back hole (see. e.g. lOhsuga et al 



re accretion inflow and outflow around the 



2009i ). that also can partly contribute to 



the width and shift of a line originated from a small disk cell. 

Concerning the disk parameters we can point out the following: (i) The 
outer radii are very similar in the sample (R ou t ~ 1500 R g ); (ii) The local 
random velocities in the disk are different from object to object and they 
are in the interval from 400 km/s to 1300 km/s. (iii) The inner radii of the 
emitting disk for all considered objects are in the interval from 300 to 600 
R 9 , (iv) The obtained inclinations are small (5° < i < 15°), and such values 
of inclinations support the idea that we are in position to more frequently 
observe the Sy 1 at face-on inclinations. 

Concerning the emission of the surrounding region we can point out that: 
i) the red-shifts are close to the cosmological, they are in the interval of ± 300 
km/s; (ii) the random velocities in this region are also different for different 
objects, and they are in the interval from 400 to 1600 km/s. 
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Figure 4: Observed averaged B./3 and Ha lines of 12 AGNs (dashed lines), fitted with the 
two-component model (solid lines) which include the emission of an accretion disk and a 
spherical emission region as a broader or double-peaked contour or a narrower Gaussian, 
respectively (dashed lines at bottom). 
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Table 1: The parameters of the disk: Zdisk is the shift and a is the Gaussian broadening 
term from disk indicating the random velocity in disk, R inn are the inner radii, R out are the 
outer radii. The zq and Wq represent the parameters of the Gaussian component. F s /Fd 
represents the ratio of the relative disk and Gaussian fluxes. The value for parameter of 
emissivity was fixed as p=3. 



Object 


% 


z disk 


a (km/s) 


Rinn \Rg) 


Rout (Rg) 


zg 


W G (km/s) 


F s /F d 


3C 120 


10 


-300 


1060 


380 


1400 


+30 


900 


0.6 


3C 273 


14 


-30 


1140 


400 


1400 


+30 


1350 


0.8 


MRK 1040 


8 


-250 


1050 


400 


1400 


00 


750 


0.2 


MRK 110 


7 


-320 


950 


500 


1400 


+ 150 


1020 


0.9 


MRK 141 


12 


-630 


1050 


300 


1100 


+200 


1620 


0.7 


MRK 493 


5 


-480 


400 


600 


1400 


+60 


360 


1.0 


MRK 817 


12 


-450 


1050 


600 


1400 


+48 


1600 


0.95 


MRK 841 


15 


-750 


1270 


450 


1400 


-300 


1500 


0.2 


NGC 3227 


12 


-780 


1080 


350 


1600 


+300 


1400 


0.8 


NGC 4253 


5 


-630 


560 


500 


1500 


-30 


600 


0.7 


PG 1116 


8 


-450 


1260 


500 


1400 


+48 


1140 


0.4 


PG 1211 


8 


-660 


830 


430 


1530 


+90 


780 


0.8 
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Table 2: The same as in Table 1, but with fixed value for parameter of emissivity as p=2.5 





% 


z disk 


/t ( l/"m / c; 1 
O ^Klll/M 


R. (R ) 

"mi \ ±L g! 


R (R } 
-Ci-out y-R-g) 


Z G 


VV g l+m/ ^ ) 


F IF, 


3C 120 


14 


+300 


960 


650 


2500 


+300 


960 


0.7 


3C 273 


18 


+90 


1191 


550 


1900 


+60 


1380 


0.9 


MRK1040 


12 





870 


700 


2700 


00 


600 


0.2 


MRK 110 


12 





890 


800 


2400 


+ 150 


1020 


1.4 


MRK 141 


14 


-630 


1040 


300 


2200 


+240 


1620 


0.8 


MRK 493 


5 


-300 


280 


800 


2400 


+60 


330 


0.9 


MRK 817 


14 


-300 


1042 


600 


2400 


+48 


1690 


1.2 


MRK 841 


23 


-300 


1270 


950 


2400 


-300 


1440 


0.3 


NGC3227 


19 


-360 


890 


800 


2600 


-300 


1500 


1.0 


NGC4253 


8 


-150 


595 


1500 


9500 


-90 


510 


0.4 


PG 1116 


15 





1480 


850 


3800 


+0 


1440 


1.0 


PG 1211 


12 


-240 


680 


730 


2530 


+90 


780 


0.9 
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Table 3: The same as in Table 1, but in search for maximal inclination 







z disk 


/t 1 [/"TY1 /CM 

O yKlLL/oj 


R- (R \ 
XL mn \ AL gJ 




Z G 


VVq ^Klll/bJ 


P 


F IF, 


3C 120 


30 


+360 


745 


1650 


19000 


+60 


960 


2.2 


0.5 


3C 273 


30 


+330 


490 


1250 


15000 


+60 


1380 


2.8 


0.9 


MRK 1040 


27 


+300 


700 


100 


17200 


00 


480 


1.3 


0.24 


MRK 110 


30 


+210 


280 


1800 


22400 


+ 150 


900 


2.0 


1.2 


MRK 141 


33 


-450 


530 


1200 


1000 


+300 


1620 


2.1 


0.8 


MRK 493 


30 


+60 


255 


10800 


124000 


+60 


270 


1.8 


0.7 


MRK 817 


35 





600 


800 


14400 


+50 


1500 


1.9 


1.0 


MRK 841 


50 


-150 


765 


1950 


27400 


-300 


1440 


2.1 


0.8 


NGC 3227 


34 


-300 


620 


1300 


12600 


-300 


1500 


2.1 


0.77 


NGC 4253 


25 


-90 


215 


2500 


69500 


-30 


510 


2.0 


0.6 


PG 1116 


30 





850 


950 


15800 


+90 


1440 


2.2 


0.3 


PG 1211 


30 





383 


920 


15530 


+90 


780 


1.9 


0.2 
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There are indications that local broadening in the disk and in surrounding 
region could be correlated (see Tables 1-3). It could indicate that those 
parameters could be linked by the same process, like wind-disk interaction, 
or something similar. 

On the other side, we compared the disk parameters obtained from the 



two-component model with the ones o btained by fitting the doub 



e-peaked 



line profiles with the disk model only. lEracleous fc Halpernl (119941 ) fitted a 
sample of 12 double peaked AGNs with the disk model. In comparison with 
their results , our disk parameters are with smaller inclination and also smaller 



outer radii. 



Strateva et al. 



(120031 ) investigated the spectra of 116 AGNs with 



double peaked lines using the disk model, and obtained inclinations of the 
accretion disks smaller than 50°, while inner radii were from 200 R 9 to 800 
R 9 and local turbulent velocities were from 780 km/s to 1800 km/s, that well 
fitted our estimated values (see Tables 1-3). Comparing those parameters 
with ones in Tables 1-3, we can conclude that they could be in a good agree- 
ment, even though our inclinations were significantly smaller than upper 
value of those given by the authors. Moreover, for double-peaked lines one 
can expect a higher inclination, because the value of inclination strongly af- 
fects the line width, even for smaller values of emissivit y index. For the outer 
radii , the estimated values were larger than 2000 R g ( lEracleous fc Halpernl . 



1994 



Strateva et al. 



20031 ). Our results sho wed smaller values o 



radii in compa r ison w ith those estimated from (lEracleous fc Halpern . 



outer 



Strateva et al. 



1994 ; 



2003J), (see Tables 1-3). 
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3. 3. Estimation of the accretion disk contribution to the BELs 

Due to a large number of free parameters, we were not able to determine 
an unique solution of parameters for the best fit (see Figs 3-6). Therefore, 
we introduced alternative approach to the analysis of the disk contribution 
to the single peak broad emission lines. 

First we set some constrains in the parameter space, and we constructed a 
grid of simulated line profiles. As a rough approximation, we fixed the width 
of the central Gaussian component to the 1000 km/s, a s mean value obtained 



from the two c ompo nent model fit (IBon et al. 



Popovic et al 



20061 ). As it was shown in 



(120041 ). the random velocities in the disk and in the non-disk 
region were approximately the same, so it was also fixed to the value of 1000 
km/s. 

The flux ratio of the simulated profiles was given by the parameter Q = 
where F s is the flux emission from the spherical region and i*disk is 



-Fdisk : 

the flux emission from the accretion disk. In the simulations the different 
values of the parameter Q have been considered, with the composite profile 
normalized to unity. The emissivity of the disk as a function of radius, r, is 
given by e = eor~ p . 

We created the grid of simulated BEL profiles with different parameters 
of the two component model. To reduce free parameters we estimated some 
constrains. Due to small change in profiles we could fix the following param- 
eters: 

• the emissivity of the disk was fixed to p=3, 

• the Gaussian broadening in the disk was fixed to 1000 km/s, 
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Figure 5: The same as in Fig. 4, but for p=2.5. 
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Figure 6: The same as in Fig. 4, but for maximal value of inclination, and with constraints 
that i > 25. 
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• the random velocities of the isotropic component was also fixed to 1000 
km/s, 

• there were no big differences in shape of the profiles if the outer radius 
was higher then several thousands R g , 

• as a mean value of results obtained from the fit, we fixed the size of 
inner and outer radius for the first sample of 12 AGN Rj„ n =400 R g and 
R O ut=3000 R g , while for the second sample of 90 SDSS AGN, these 
parameters were fixed by averaging the disk si ze obtained from the 



fitting of the double-peaked lines in the work of lEracleous fc Halpern 



(120031 ) . so we fixed the inner radius to be Rj nn =600 R g , and the outer 



R out =4000 R g . 



It was noticed that for the change of the shape of BEL profiles, variation 
of inclination and the flux ratio played dominant role, so the grid of models 
was finally reduced only to these two parameters. Simulated spectra were 
compared to the spectra from different samples. 

First the profiles were normalized to the same intensity. We noticed that 

)e neglected since normalize d 



Bon et al. 



2009 



Bon . 



20081 ). 



the variation of the emissivity parameter could 
profiles showed very small change in shape (see 
For that reason, this parameter was fixed to the value p=3. 

For each line profile, the full widths at 10%, 20%, 30% of maximum in- 
tensity were measured, and compared after being normalized to full width at 
half maximum, which was introduced through coefficients kj (j = 10, 20, 30) 
as k w = wi %/u>50%, k 20 = w 20 %/w 50 % and k 30 = w 30 %/w 50 %. 

Fixing the inner and outer radius to an averaged value, obtained from 
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the study of BELs with double-peaked lines (lEracleous fc Halpernl . 120031 ) . 
we estimated the values of i and Q for all analyzed galaxies. 

For the first sample, the inner radius of the disk was chosen to be Rj nn — 



400 R n, as the mean value of this parameter in the previous fittings ((Bon 



20081 ). For the outer radius, the chosen value was R, 



out 



3000 R G , al- 



though the mean value of the outer radius for the sample was about 30000 
( )Bon et all 120061 ). since the variation of this parameter, after values of sev- 
eral thousand gravitational radii, did not show significant influence to the 
shape of the line (especially if the line is noisy). 



4. Results 

In order to estimate possible disk presence in the single peaked broad 
emission lines of AGN, we applied several methods. 

From the Gaussian analysis we could conclude that two distinct kinematic 
regions might be present in the BELs. The central Gaussian component 
indicated the existence of the emission from the region with intermediate 
velocities (ILR), while in the wings two shifted broad Gaussians indicated 
possible presence of disk-like emission (VBLR). The results for the widths of 
the emission lines of Fe II template showed that these lines probably originate 
in the ILR (see Figures 1-3). 

From the fit with two component model where one component represents 
the emission from the accretion disk (mostly affects the wings of the emis- 
sion line), and another from the surrounding region with isotropic velocity 
distribution of the emitters (that contribute to the core of the line), we can 
conclude: 
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Figure 7: The measured width ratios (crosses-sample of 90 SDSS, asterisks-sample of 12 
H a , x signs-sample of 12 Hp) and simulated values (dashed lines) from the two-component 
model for the different contribution of the disk emission to the total line flux (Q=0.3, 0.5, 
0.7, 1, 1.5 and 2) and different inclinations are considered (solid iso lines presented i = 
10, 15, 20, 25 and 30 degrees, respectively) /JQ 



• two component model can well fit the profiles of broad emission lines 
of our samples, but it is very hard to determine the disk parameters, 
because of the large number of the free parameters (see Figs 4-6), 

• random velocities in the surrounding region with isotropic velocities are 
similar to random velocities in the disk, what implies that these two 
regions could be connected through some process (for example through 
the wind of the accretion disk), 

• the values of the parameters indicate that the inclinations are smaller 
than 50°, and the dimensions of the disk could be from several hundred 
R g for the inner radius, to several hundred thousands R g for the outer 
radius. 



After comparing the measurements of ki for the sample, with the mea- 
surements of the simulated profiles, we could derive the estimates of i and 
Q. As the result, we found that the most of the measured points were be- 
tween 0.5 < Q < 1.5, and 10° < i < 25° (see Fig. [7]). The results showed a 
consistence between estimations of i and Q fo r both Ha and H/?, with very 



small discrepancies (see Table 1 in 



Bon 



2008|). 



After this analysis , we tested this method to the larger sample of 90 AGN 



sec 



Bon et al. 



2009i ). The measurements were performed in the same way 



as with 12 AGN, with the difference that we compared it to the simulated 
spectra with the disk dimensions of Ri nn = 600 Rg and R out = 4000 Rg 
(averaging the disk size obtained from the fitting of the double-peaked lines 
in the work of (jEracleous fc Halpernl . 120031 )). With this sample we measured 
only Ha lines. As it could be seen in Fig. [71 most of the measurements are 
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located within 1 < Q < 2, and 10° < % < 25° (see 



Bon et al 



2009). 



5. Conclusions 

In this paper we outline our recent investigation about the presence of 
the disk emission in single peaked BELs. 

We carried out the following analyzes: i) fitting of single peaked BELs 
with Gaussians and with two component model, and ii) comparing line pa- 
rameters of a sample with simulated profiles. 

From our investigations we can conclude that there could be a high prob- 
ability that the disk emission flux might be present in the single peaked 
emission line profiles. In principle, the contribution of the disk emission to 
the total flux of the single peaked BELs is in most cases smaller then 50%. 
Also, we found that the disk inclinations were mainly smaller then i < 25°. 
Such small inclinations and large discrepancies from simulated profiles for 
higher inclinations should be discussed in the contest of unified model, i.e. 
possible disk orientation to the torus or partial obscuration by the torus. 
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